Partial dissociation of water on Fe 3 O4(001): adsorbate induced charge and orbital 
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The interaction of water with Fe3O4(001) is studied by density functional theory (DFT) calcula- 
tions including an on-site Coulomb term. For isolated molecules dissociative adsorption is strongly 
promoted at surface defect sites, while at higher coverages a hydrogen-bonded network forms with 
alternating molecular and dissociated species. This mixed adsorption mode and a suppression of 
the (y/2 x y/2)R45° -reconstruction are confirmed by a quantitative low energy electron diffraction 
(LEED) analysis. Adsorbate induced electron transfer processes add a new dimension towards 
understanding the catalytic activity of magnetite(OOl). 
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Magnetite has attracted continued interest in the past 
decades due to its fascinating properties: The high tem- 
peraturephase is predicted to be a half-metallic ferri- 
magnet [l[ with high magnetic ordering temperature. At 
around 120 K it undergoes the so called Verwey transi- 
tion ■ Both the type of transition (metal-to-insulator 
vs. semiconductor-to-semiconductor [H, HJ) and the type 
of charge (CO) and orbital ordering (00) at the octahe- 
dral iron sites in the low temperature phase are subject 
of an ongoing debate [j| H, 0, H] • 

Besides its applications in magnetic recording and as a 
prospective material for spintronics, magnetite acts as a 
catalyst e.g. in environmental redox reactions [i| [l(| , or 
in the water gas phase shift reaction [ll[ . Typically these 
reactions take place in an aqueous environment prompt- 
ing the need to understand how water interacts with the 
Fe304-surface. Water can bind to a surface in different 
modes (e.g. molecular or dissociative) leading to a vari- 
ety of functional groups that can affect significantly the 
surface properties and availability of reaction sites and 
result in a complex surface chemistry fl2l . Il3j . 

The catalytic activity of magnetite is typically re- 
lated to the presence of both ferrous and ferric iron in 
its inverse spinel structure. In the [001]-direction two 
types of layers alternate: A-layers with tetrahedral iron 
(FeA 3+ ) and B-layers containing oxygen and octahedral 
iron (Fc B 2+ ' 3+ ). Both bulk truncations of Fe 3 O 4 (001), 
either with an A- or a B-layer are polar of type three 
according to the classification of Tasker [l4| . To explain 
the origin of the (V2 x v / 2)i?45° -reconstruction previ- 
ous surface models proposed an ordering of surface de- 
fects (e.g. (TSUI)- Recently, DFT calculations [LI G3 
have shown that the symmetry lowering at Fe3O4(001) 
is achieved rather by a distorted B-layer, supported also 
by surface x-ray diffraction [17|. LEED [19j and scanning 
tunneling microscopy (STM) [201 ] . 
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Despite its importance, only a few studies have ad- 
dressed the interaction of water with Fe3O4(001). While 
initial adsorption was related to surface defect sites, an 
extensive hydroxylation was reported from x-ray photoe- 
mission experiments (XPS) beyond a threshold pressure 
of 10~ 3 — 10~ 4 mbar[2lj. Both this study and temper- 
ature programmed desorption (TPD) measurements [12] 
indicate multiple adsorbate sites on the surface. Molec- 
ular dynamics simulations with empirical potentials [23], 
[2~ij point towards a dissociative adsorption. However, 
such studies cannot provide reliable information on the 
energetics and underlying electronic phenomena. 

Here we address these fundamental questions within 
a combined DFT and LEED approach. By varying the 
coverage and adsorbate configuration of water molecules 
we compile a surface phase diagram in the framework 
of ab initio thermodynamics. The results show that iso- 
lated molecules dissociate on the clean Fe3O 4 (001) sur- 
face. This process is strongly favored at oxygen vacan- 
cies. With increasing coverage a crossover to a mixed 
adsorption (both molecular and dissociative) takes place 
that is confirmed in a quantitative LEED analysis. LEED 
shows also a strong suppression of the (y/2 x y/2)m5° 
surface reconstruction. Furthermore we find that the ad- 
sorbed species invoke electron transfer processes in the 
subsurface layers resulting in an unique charge and or- 
bital order that may have important implications on the 
catalytic activity of the surface. 

DFT calculations were performed using the full po- 
tential linear augmented plane wave method in the 
WIEN2K [25( implementation. The generalised gradi- 
ent approximation (GGA-PBE) [26| of the exchange- 
correlation potential is used. Systematic studies [U EH 
show that the PBE functional gives a good description of 
hydrogen bonding w.r.t. equilibrium geometries. The er- 
ror of overbinding which can be as large as 20 meV/bond 
will not affect our conclusions. Because magnetite is a 
strongly correlated material we have explored the influ- 
ence of electronic correlations within LDA/GGA+J7 [2^| 
using U = 5 eV and J = 1 eV [3(| , similar to values used 
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for bulk Fe3C>4 0, 0, [H|. The surfaces were modelled 
in the supercell geometry with slabs containing seven B- 
layers and six A-layers separated by a vacuum region of 
10-12 A[33]. The positions of adsorbates and the atoms 
in the outer two BA-layers were relaxed. 

We have compared the stability of more than 
30 different configurations using ab initio ther- 
modynamics [33|]. The surface energy, 7(T,p) = 

TK (^Fe O 3 O 4 (00i) ~ ^FcMFc — ^o^o ~ ^Vh 2 o^h 2 o) de- 
pends on the Gibbs free energy and the chemical 
potentials of the constituents. The Gibbs free energy 
can be expressed through the total energy from the 
DFT-calculations (33|. Because there are two species in 
the gas phase, O2 and H2O, the surface phase diagram 
is three dimensional. Fig. [T] displays a two dimensional 
projection with the most stable configurations for given 
/^o, A i H 2 o- We first consider the termination of the 
clean surface as a function of oxygen pressure. As 
discussed previously, a modified B-layer (denoted as B) 
showing lateral and vertical distortions in the surface 
layer with a wave-like pattern is favored over a broad 
range of oxygen pressures. However, at oxygen poor 
conditions a B-layer with oxygen vacancies (B+Vo) is 
stabilized. This defective surface, previously proposed 
in a STM study [Hj], shows dramatic relaxations where 
the oxygen opposite the vacancy moves towards the 
Fee-row by ^0.8 A. Starting from these two termina- 
tions [34J, we have adsorbed water molecules, varying 
their concentration and geometry. We find a strong 
tendency for isolated molecules to dissociate in surface 
oxygen vacancies (Vo+OH) whereby a proton diffuses 
to a surface oxygen (0 s ) further away and forms a 
surface OH-group (O s +H). Thus, even at very low water 
vapor pressures, all surface defects are expected to be 
filled with OH-groups, consistent with XPS results [2lj . 
The dissociation of water in the defect sites invokes 
significant changes in the electronic structure as shown 
in Fig. [2k and c, involving the In and 3cr molecular 
orbitals of the OH groups and a switching from Fe 3+ to 
Fe 2+ of Fe B S-l underneath O s +H. 

The adsorbate-adsorbate interaction that sets in, 
when a second molecule is adsorbed within the (\/2 x 
%/2)-R45° unit cell, leads to a mixed adsorption mode: 
One molecule dissociates protonating a surface oxygen, 
while the intact molecule forms a hydrogen bond with 
the OH-group. The main part of the phase diagram is 
dominated by a mixed adsorption mode of four water 
molecules where all surface Fes-sites (FeeS) are satu- 
rated (B+4H 2 0(M)). Full dissociation is 12 meV/A 2 
less stable. Full hydroxylation of the surface is also not 
likely to occur as the formation of a surface OH-group 
neighboring a subsurface FeA is extremely unfavorable. 

The electron density redistribution with respect to the 
B-layer in Fig. [2p indicates a weak charge accumula- 
tion between FeeS and the water molecule (O w ), but the 
strongest charge rearrangement takes place between Fee 
and the OH-group (0 OH ) of the dissociated molecule, in- 
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FIG. 1: Bottom view of the surface phase diagram of 
Fe3O4(001) showing the most stable configurations for given 
(/io,MH 2 o): B+Vo (red), B-layer (magenta), dissociated 
molecule in an oxygen vacancy B+Vo + 1H20(D) (grey 
hatched semitransparent area), mixed adsorption of two 
(green) and four water molecules on a B-layer (purple). The 
ranges of and /j,h 2 o correspond to the vapor phase of H2O. 
\io and /xh 2 o have been converted into pressures for 300 and 
800 K. Additionally, the top views of the most stable con- 
figurations are displayed with positions of O, FeB, FeA and 
H marked by cyan, purple, orange and white circles, respec- 
tively. In B+Vo a white cross marks the position of the va- 
cancy. 



volving depletion of 3ct and accumulation in 1-7T molecular 
orbitals of OH along with depletion of the d z 2 orbital at 
FeeS. The main actuator of partial dissociation appears 
to be the formation of strong intermolecular hydrogen 
bonds. As a result both H2O and OH tilt from the on- 
top position towards each other and O w -H elongates from 
0.95 A (gas phase) to 1.10 A, resulting in an O w -H- • -0 OH 
of 2.47 A. A fingerprint of the mixed adsorption are the 
two distinct bond lengths between Fee-O w and Fee-0 
of 2.09 A and 1.97 A respectively. This feature is con- 
firmed by the LEED analysis described below. 

LEED measurements were performed on a synthetic 
magnetite sample. The clean surface, prepared by Ar + - 
ion sputtering and subsequent annealing at 900-1000 
K at po 2 = 5 x 10~ 7 mbar for 2-3 h, exhibited a 
(\/2 x v / 2)^?45°-LEED pattern with sharp superstruc- 
ture reflections and low background (Fig. [3] a). Water 
was adsorbed at 273 K for 2 min at a water vapor pres- 
sure of 2 x 10 -6 mbar, resulting in a (1 x 1)-LEED pattern 
with an enhanced background (Fig. [3Jd). While anneal- 
ing up to 770 K restores the superstructure spots, the 
shape of the the LEED I(V) curves cannot be recov- 
ered 37], indicating some residual hydroxylation of the 
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FIG. 2: Adsorbate induced electron density redistribution for 
(a) B+Vo+lH 2 0(D) and (b) B+4H 2 0(M). Red/blue corre- 
sponds to regions of charge accummulation/depletion. Rele- 
vant bond lengths are given in A. Total and projected density 
of states (DOS) before (dashed line) and after water adsorp- 
tion (solid line, grey area) for (c) B+Vo + lH20(D) and (d) 
B+4H2 0(M). S (S-l) denote surface (subsurface) Fee sites. 



surface, as observed by Kendelewicz et al. [21| . Seven 
LEED I(V) curves of the (1 x l)-surface were measured 
in the energy range 50 - 300 eV at 100 K. An increase 
of the background current during the measurement is at- 
tributed to uncorrelated defects induced by the electron 
beam 38]. However, no changes were observed in the 
LEED I(V) curves during repeated measurements over 
periods of several hours. 

LEED calculations were performed with the layer dou- 
bling method and a least squares optimization [35| us- 
ing constraints for bond lengths. The crystal potential 
was calculated from a superposition of atomic potentials 
using optimized muffin-tin radii [36| which led to reli- 
able structural determination of the clean FesO^OOl)- 
surface 10 phase shifts were used. All positions and 
occupation numbers within the adsorbate and top B-A-B 
substrate layers were optimized in a (V2x v / 2)i?45° -unit 
cell, while thermal parameters were kept fixed. The best 
fit (Rp = 0.27) was obtained for a model where all FeeS 
sites have adsorbed oxygen on top. Both the surface 
Fee- and adsorbate sites are occupied by ~ 80% possibly 
due to defect creation during the preparation procedure. 
The adsorbed O shows strong lateral shifts by ~ 0.4 A off 
the Fee sites in agreement with the DFT results (0.23- 
0.28 A). The main feature are the two different bond 
lengths Fe B -O w /0 OH : 2.12 A and 1.93 A, confirming the 
simultaneous occurrence of hydroxyl groups and molecu- 
lar adsorption. Further details on the structural analysis 
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FIG. 3: (Color online) LEED pattern a) before and b) after 
water adsorption. The superstructure spots in a) indicative 
of the (y/2 x \/2)^45 surface reconstruction are largely sup- 
pressed in b). c) experimental (averaged over symmetrically 
equivalent beams - black) and calculated (red/grey) LEED 
I(V) curves of the water adsorbed surface. 



will be published elsewhere |37j |. 

We turn next to the surface and adsorbate induced 
electronic effects on Fe3O4(001). The electron density 
plot in Fig. 3] displays the minority ti g occupancies at 
the Fee-sites, thus allowing to distinguish between sites 
with predominant Fe 2+ or Fe 3+ -character [3(J. The no- 
tation Fe 2+ and Fe 3+ is used here for simplicity, but 
the difference in total 3d occupation within the muffin- 
tin sphere is much smaller (0.2-0.4e) consistent with 
XRD- [5[ and LDA+C7 studies for the low tempera- 
ture bulk phase @, Q. The magnetic moments al- 
low a clearer discrimination: M Fc 2+ = 3.54 — 3.75/xb, 
M Fe 3+ = 3.90 — 4.10/ib, respectively. For most of the 
stable systems we find that the surface layer contains 
exclusively Fe 3+ while a unique charge and orbitally or- 
dered state emerges in the deeper layers depending on 
the type of termination In B+V (Fig. Hi) the 

two subsurface Fee next to Vo are Fe 2+ (d xz ± d yz ). 
Upon water adsorption the positions of Fe 2+ switch in 
B+Vo + lH20(D): now the two Fee below the surface 
OH groups are Fe 2+ , resulting in alternating Fe 2+ , Fe 3+ - 
sites. In the mixed adsorption case (B-r-4H 2 0(M)) the 
charge ordering in the subsurface layer is not significantly 
altered compared to B-layer, where one out of four FesS- 
1 is Fe 2+ and a second is in an intermediate valence 
state. However, the slight tilting of the t2 g orbital at 
Fee 2+ induces a completely different OO in S-2: from 
d xz ± d yz to alternating d xy and cL^-orbitals. Recent 
GGA+U calculations have shown that in bulk magnetite 
a variety of CO/OO states can be realized by symme- 
try lowering j31(. Here the presence of a surface and 
adsorbates imposes a unique CO/OO state for each ter- 
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FIG. 4: Side view of the electron density integrated be- 
tween -1.3 eV and Ef showing the occupation of the minor- 
ity t2 g orbitals at the FeB-sites. a) B+Vo; b) B-layer; c) 
B+Vo+lH 2 0(D); d) B+4H 2 0(M). For the colour code see 
Fig.rj] 



mination reaching several layers below the surface. More- 
over, most of them violate the Anderson criterion: Fe404 



cubes in the surface and subsurface layer are predomi- 
nantly electron poor (3Fe 3+ +lFe 2+ ), while those in S-l 
and S-2-layer are electron rich (lFe 3+ +3Fe 2+ ). Evidence 
for orbital ordering in thin magnetite films was reported 
recently from resonant soft x-ray diffraction [8j]. 

The charge ordered state in the B-layer leads to the 
opening of a band gap of 0.3 eV as shown in Fig. 
consistent with previous calculations II 811 and scanning 
tunneling spectroscopy measurements [40( . However the 
adsorption of water on the surface and the formation of 
surface hydroxyl groups leads to half-metallic behavior. 

We have shown that on FesO^OOl) water tends to dis- 
sociate in oxygen vacancies or partially dissociate on the 
nondefective surface in contrast to the full dissociation 
found on FesO^lll) [2l|, [4l[. This adsorption mode is 
triggered by the presence of both Lewis acid and base 
sites on the surface and the intermolecular interaction as 
reported also on rutile(llO) [i^ ]. anatase(OOl) [43| and 
MgO(001)-surfaces Q. The hydrogen bonded OH and 
H2O can easily rearrange and thus provide adsorption 
sites for further species (e.g. heavy metal complexes). 
Our results indicate a pathway to manipulate e.g. the 
catalytic properties of transition metal oxide surfaces by 
triggering electron transfer processes and inducing new 
charge and orbitally ordered states via adsorbates. 
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